arXiv: astro-ph/0112410v2 20 Dec 2001 


Mon. Not. R. Astron. Soc. 000 , 000-000 (0000) 


Printed 1 February 2008 (MN LMjrjX style file v2.2) 


Shocked molecular gas towards the SNR G359.1—0.5 and 
the Snake 

J. S. Lazendic, 1,2 * * M. Wardle, 1 M. G. Burton , 3 F. Yusef-Zadeh , 4 J. B. Whiteoak , 2 
A. J. Green , 1 M. C. B. Ashley 3 

1 School of Physics A28, University of Sydney, Sydney NSW 2006, Australia 

2 Australia Telescope National Facility, CSIRO, PO Box 76, Epping NSW 1710, Australia 

3 School of Physics, University of New South Wales, Sydney NSW 2052, Australia 

4 Department of Physics and Astronomy , Northwestern University, Dearborn Observatory, 2131 North Sheridan Road, Evanston, 

IL 60201-2900 


Accepted 


ABSTRACT 

We have found a bar of shocked molecular hydrogen (H 2 ) towards the OH(1720 MHz) 
maser located at the projected intersection of supernova remnant (SNR) G359.1-0.5 
and the nonthermal radio filament, known as the Snake. The H 2 bar is well aligned 
with the SNR shell and almost perpendicular to the Snake. The OH(1720 MHz) maser 
is located inside the sharp western edge of the H 2 emission, which is consistent with 
the scenario in which the SNR drives a shock into a molecular cloud at that location. 
The spectral line profiles of 12 CO, HCO + and CS towards the maser show broad-line 
absorption, which is absent in the 13 CO spectra and most probably originates from 
the pre-shock gas. A density gradient is present across the region and is consistent 
with the passage of the SNR shock while the H 2 filament is located at the boundary 
between the pre-shocked and post-shock regions. 

Key words: supernova remnants - ISM: clouds, masers, shock waves - individual: 
G359.1-0.5, Snake - infrared: ISM 


1 INTRODUCTION 

OH(1720 MHz) maser emission detected towards supernova 
remnants (SNRs) has been attri buted to shock waves drive n 
into adjacent molecular clouds ( Frail, Goss fe Slysh 1994 ). 
The production of this maser emission, without the simul¬ 
taneous production of maser emission in the other three 


along the radio-continuum sh ell of the SNR G359.1-0.5 
(Yusef-Zadeh et al. 1995, [19961 ). The b rightest of the masers 
(maser A from Yusef-Zadeh et al. 1995) occurs where a non¬ 
thermal filament known as ‘the Snake’ appears to cross the 
western edge of the SNR shell. The filament differs from 
other Galactic Centre nonthermal filaments i n showing a 


ground istafu OH transitions at 1012, 1005 and 1007 MHz, 
require: specific cuudiliuus (gas density n\i 2 < MCf' uiU'\ gas 

temper ature Tfc m ^50 - 125 K, OH column density Auh 1 u 

10 15 - 10 16 cm -2 ) which can be attain ed in the cooling gas 
behind a non-dissociative shock wave ( Lockett J _Gauthier& 
Elitzui_JT99 Wardle 1999 ). Surveys by Frail et al. (1996) 


number of kinks along its 20 aremin extent ( Gray_eyjd 


^L99f ), and it s origin is not well established ( Nicholls & Le 
Strange lO gSl^^eniordMOTZ^Uchida et al. 1996). The latest 


Green et al. (1991) and [Koralesky et al. (1998) have found 


that about 10 pci cent uf 'M80 ubscivcd Galactic 3NRs 
are ass ociated with OH(1720 MHz) mascra. The detection 
fraction increa ses to 30 per cent for obje cts in the Galactic 
Centre region (Yusef-Zadeh et al. 1996a). 

Perhaps the most striking result of the survey towards 
the Galactic Centre region is the detection of OH masers 


* Current address: Harvard—Smithsonian Center for Astro¬ 
physics, 60 Garden Street, Cambridge, MA 02138, USA 


model by p3icknell fc Li (2001 ) proposes that the Snake is a 
magnetic flux tube anchored in dense rotating material. It 
has been suggested that the Snake is interacting with the 
shell of G359.1-0.5 because both objects show change in 
brightness at the apparent crossing point (Uchida, Morris 


& Yusef-Zadeh 1992; Gray et al. 1995). This could, how¬ 


ever, be an artifact due to superposition of disparate com¬ 
ponents along the line of sight. Nevertheless, the presence of 
an OH(1720 MHz) maser - a signature of SNR/molecular 
cloud interaction - at this location is quite intriguing. In 
order to investigate the proposed interaction between the 
Snake and the SNR, and to characterise the ambient molec¬ 
ular gas in the region, we observed molecular hydrogen (H 2 ) 
and a number of other molecular species, particularly search- 
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ing for signatures of shocked gas. The observations are dis¬ 
cussed in section 2, the results are presented in section 3 
and discussed in section 4. The conclusions are presented in 
section 5. 


2 OBSERVATIONS 


2.1 UNSWIRF observations 

The observations towards maser A (RA(1950) = 

17 h 41 m 46 s 042, Dec. (1950) = -29° 49' 51'.'03) were carried 
out in June 1998 with UNSWIRF^] on the Anglo-Australian 
Telescope (AAT). Two data sets comprising five and six 
Fabry-Perot frames, each spaced by 40kms _1 , were taken 
in the 2.12 fin i H 2 1-0 S(l) transition. An additional frame 
was taken at a setting of —400kins -1 to enable subsequent 
continuum subtraction. The integration time was 120 sec¬ 
onds per frame. A resulting image has diameter of 100 arcsec 
and a pixel size of 0.77 arcsec. The velocity resolution was ~ 
75 kins -1 FWHM. All the data were reduced using m odified 


routines in t he iraf software package as described in Ryder 
et al. (1998). Intensity calibration was performed using the 
standard star BS 5699. A velocity cube was constructed by 
averaging the frames at the same Fabry-Perot setting from 
two data sets, which was then fitted with the instrumental 
Lorentzian profile to determine the line flux (to within 30 
per cent) and central velocity (to within 20kms~ 1 ) for the 
H 2 emission across the field. We note that these data do not 
provide a true velocity-spatial cube of the emission, as the 
instrument does not have the velocity resolution to provide 
an accurate line velocity at each position. Rather, the line 
velocity is the central velocity at each position. Changes in 
line center velocity can, however, be determined. To con¬ 
firm the detection, the observations were repeated in June 
1999 for the 1-0 S(l) line and additional observations were 
obtained in the 2.25/nn H 2 2-1 S(l) line. On this occasion 
five Fabry-Perot frames were used for both transitions, with 
an integration time of 180 seconds per frame. To establish 
the coordinate scale for the UNSWIRF images, we used the 
positions of stellar sources in the K-band continuum image 
listed in the Two Micron All Sky Survey (2MASS) point 
source catalogue (Cutri 1997) (which has an accuracy of 
«0"l). 


2.2 SEST observations 

A millimetre-line survey towards maser A was carried out 
with the 15-m SEST^j at La Silla, Chile during June 2000. 
Position-switching was used with a reference position at 
RA(1950) = 17 h 45 m 00 s , Dec. (1950) = -31° 00' 00". For 
periodic pointing calibration and focusing the 3-mm SiO 
maser of AH Sco was observed. The main-beam brightness 


1 Unive rsity of New South Wales Infrared Fabry-Perot (Ryder et 
al. 1998| ) 

- The Swedish-ESO Submillimetre Telescope (SEST) is operated 


by the Swedish National Facility for Radio Astronomy, Onsala 
Space Observatory and by the European Southern Observatory 
(ESO). 


Table 1. The first column lists observed molecular species, the 
next three columns list the SEST observational parameters and 
the last column gives observed main—beam brightness tempera¬ 
ture (T m j,) and rms noise towards the maser A. 


Molecule 

Transition 

Frequency 

(GHz) 

Beam Size 
(arcsec) 

Tmb 

(K) 

13 co 

2-1 

220.399 

23 

0.78T0.21 


1-0 

110.201 

45 

1.20F0.26 

12 co 

2-1 

230.538 

23 

16.00i0.20 


1-0 

115.271 

45 

10.30i0.50 

CS 

3-2 

146.969 

34 

0.42i0.12 


2-1 

97.981 

52 

0.54i0.14 

HCO+ 

1-0 

89.188 

54 

0.48i0.08 

O 

00 

0 

2-1 

219.560 

23 

< 0.16 


1-0 

109.782 

45 

< 0.12 

HCN 

1-0 

88.632 

55 

<0.17 

H 2 CO 

3(2,2) _ 2( 2i l) 

218.475 

24 

< 0.13 


3(0,3) — 2(0,2) 

218.222 

24 

< 0.15 

SiO 

5-4 

217.106 

24 

< 0.14 


2-1 

86.848 

57 

<0.17 


The first part of the table lists results from Gaussian fit of the 
two 13 CO transitions, whose line profiles peak at — 13kms _1 and 
have line widths of ~15 - 20 km s 1 . The second part of the table 
lists the parameters of 12 CO, CS and HCO + transitions. Their 
line profiles suffer broad-line absorption and Gaussian fitting was 
not possible. We therefore list only the maximum T m /, of their line 
profiles, which occur at — 5kms — T The third part of the table 
lists the non—detections. Upper limits are 2a values. 


temperature (T m f,) scale of the spectra was calibrated on¬ 
line using a black-body calibration source at ambient tem¬ 
perature, and later corrected for the telescope main-beam 
efficiency (0.74, 0.70, 0.67 and 0.45 at 85 - 100 GHz, 100 - 
115 GHz, 130 - 150 GHz and 220 - 265 GHz respectively). Si¬ 
multaneous observation of transitions in two different bands 
was available using dual SIS receivers at 3 and 2 mm, or 3 
and 1.3 mm. 

For spectral- line observing we used an acousto-optical 
spectrometer split into two 43-MHz bands, each with 1000 
channels, providing velocity coverages ranging from 60 to 
100 km s -1 and velocity resolution ranging from 0.06 to 
0.14kms^ 1 across the total wavelength range. A low- 
resolution spectrometer with resolution of ~l-2kms~ 1 was 
also used simultaneously to determine the spectral base¬ 
lines in the case of broad-line emission. The final data were 
smoothed over three channels. The observed transitions, fre¬ 
quencies and antenna beam sizes (FWHM) are listed in Ta¬ 
ble |l]. For the 12 CO, 13 CO and CS transitions we obtained 
maps with approximately half-beam and one-beam sam¬ 
pling intervals at lower and higher transitions respectively. 
The final images covered the 2x2 arcmin 2 region around 
maser A that coincides with the UNSWIRF field of view. 
The integration time was 30 seconds for CO isotopomers 
and 90 seconds for CS transitions. A set of spectra on-source 
and at offsets of 20 arcsec in RA and Dec. was obtained for 
each of the other molecular transitions listed in Table [jj with 
integration times of 180 seconds at each position. 
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Figure 1. Image of the velocity-integrated 2.12 pm H 2 1 
0 S(l) line (contours) overlaid onto the 2MASS K-band im¬ 
age (greyscale). The contours are: 13, 20, 30, 45, 56 and 
65 X10 -6 erg s —1 cm -2 sr~ 1 . The cross marks the location of 
OH(1720 MHz) maser A. 


Figure 2. The line-centre velocity distribution of the 2.12 /im 
H 2 1—0 S(l) line emission in both greyscale and contours. 
The greyscale units are in kms' 1 and the contours are: 
—5, —10, —15, —25, —35 and —55kms —1 . The cross marks the lo¬ 
cation of OH(1720 MHz) maser A. 
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3 RESULTS AND ANALYSIS 
3.1 H 2 emission 

Fig. |l] shows contours of velocity-integrated 2.12 (im H 2 
1-0 S(l) line emission superimposed on a greyscale image 
of the 2MASS K-band. Some small-scale structure appar¬ 
ent in the H 2 distribution are residuals from subtracting 
the stellar continuum. A bar of H 2 emission, 1.5 arcmin 
in length and 15 arcsec in width, extends from the north¬ 
east to the south-west. We note, however, that there is 
weak extensive diffuse emission across the observed region, 
which was not included in fitting procedure because low 
S/N. The lowest contour in Fig. M thus corresponds to a 
(4 a) flux level of 2x 10 -5 erg s -1 cnh 2 sr -1 . The peak flux 
density is 10.6xl0 -4 ergs -1 cm -2 sr -1 and the total flux 
density is 9.2xl0 -13 ergs -1 cm -2 . After correcting for the 
typical extinction of the Galactic Centre region (3 mag¬ 
nitudes in K-band) the parameters of the 1-0 emission 
are: peak flux density of 1.7xl0 -2 ergs -1 cm' 2 sr” 1 , total 
flux density of (1.4±0.4) x 10 -11 erg s -1 cm -2 and luminos¬ 
ity of ~33 Lq. There is no obvious emission in the indi¬ 
vidual frames of the 2.25 pm H 2 2-1 S(l) data. However, 
adding the pixels in the central frame of the 2-1 cube over 
the region delineated by the 1-0 emission, we derive a to¬ 
tal flux density of (2.4±0.4) x 10 -14 ergs -1 cm -2 , which is 
~ 2.9x 10 -13 ergs -1 cm -2 after correcting for extinction. 
The corresponding total flux density in the 1-0 peak frame 
is 3.1xl0 -13 ergs -1 cm -2 , i.e. 4.9xl0 -12 ergs -1 cm -2 after 
correcting for extinction. The ratio between the 1-0 and 2-1 
line emission is thus ~ 20. 

The line-centre velocity distribution of H 2 1-0 S(l) 



Figure 3. A cut, 5 arcsec wide, along the centre of the H 2 bar 
extending 80 arcsec diagonally from the top of the bar. 


shows a cut along the H 2 bar illustrating the velocity gradi¬ 
ent. Note that the small-scale structure showing deviations 
of ~10-15 kms -1 probably reflects the errors in determin¬ 
ing the velocity at each position. 

The presence of a star (RA(1950) = 17 h 41 m 46)5 , 
Dec.(1950) = —29° 49' 45V0) in the centre of the H 2 fil¬ 
ament, together with the velocity difference between the 
north-east and south-west sections, and the jet-like appear¬ 
ance of the brightest H 2 contours would be consistent with 
a collimated bipolar jet originating from the star. However, 
other considerations suggest that this is unlikely. Firstly, 
the star in question cannot be a protostar. The 2MASS 
point source catalogue gives the J, H and K magnitudes 
of the star as 13.22 ± 0.04, 10.31 ± 0.05 and 8.97 ± 0.04 re¬ 


spectively. Using the JHK intrinsic colours from Bessell & 


comparable with the average values of about —13 kms -1 for 
the northern part and —35kms -1 for the southern part ob¬ 
tained from the repeated observations in June 2000. Fig. N 


the H 2 rar, with the mean velocity ranging from Okms 1 in 

Brett (198$) and 

Koornneef (1983 

), and the reddening vec- 

the north to —50kms -1 in the south. These velocities are 

tor from 

Bessell & Brett (1988), we find that the extinction 


in K-band is Ak ~ 2 (visual extinction Av ~ 20). These 
values are consistent with the lack of detection in the Digi¬ 
tized Sky Surveys (DSS). After correcting for extinction, we 
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Figure 4. The velocty-integrated image of the H 2 line emis¬ 
sion (contours) and the 6-cm radio continuum image (greyscale). 
The H 2 filament is located where the Snake (the dark east-west 
bar) crosses the shell of G359.1-0.5 (dark band extending almost 
north-south). The greyscale units are in Jy beam - 1 and the con¬ 
tours are: 13, 20, 30, 45, 56 and 65 xlO -6 ergs -1 cm -2 sr -1 . The 
cross marks the location of OH(1720 MHz) maser A. 


find that derived magnitudes are consistent with either a K 
dwarf within 100 pc of the Sun, or a K or M giant at a dis¬ 
tance of several kpc. The former case can be excluded given 
that there are no molecular clouds within 200 pc of the Sun 
that could give rise to the H 2 emission. In the latter case 
the bright star is evolved and cannot be driving the outflow. 
The apparent location of the star at the centre of the H 2 bar 
is, therefore, most probably a chance alignment. Of course, 
it is possible that apparent outflow is being produced by a 
star at several kpc which is too faint to be detected. A more 
serious problem is that the emission does not resemble that 
of other outflows detected in the 1-0 S(l) line . In bipolar 
outflows, such as Cep E (Eisloffel et al. 1996), the emis¬ 
sion associated with jets generally peaks in discrete knots 
or lobes at some distance from the driving source. Further¬ 
more, the jets are narrow close to the driving source and the 
lobes have distinct, roughly constant velocities. The veloc¬ 
ity gradient and morphology, including the extended diffuse 
component, present in our H 2 data are thus not consistent 
with entrainment by a fast, narrow protostellar jets. 

Fig. y shows contours of velocity-integrated H 2 line 
emission superimposed on a 6-cm greyscale VLA image with 
a resolution of 12 "7 x 8 "2 (P.A. = 62°) ( Yusef-Zadeh et 
al. 200||). The bar of H 2 emission is almost perpendicular 
to the Snake and it runs parallel to the edge of an elliptical 
depression located ~15 arcsec to the west of the H 2 bar. 
The synchrotron emissivity of this depression is reduced by 
< 30 per cent and it surrounds the Snake filament where it 
crosses the shell of G359.1-0.5. 


3.2 Characteristics of associated molecular cloud 

3.2.1 Identifying the cloud associated with G359.1-0.5 

Observations of molecular clouds associated with SNRs are 
potentially subject to confusion due to other molecular 


clouds along the line of sight. However, the maser veloc¬ 
ity of -5kms _1 can be related to the systemic velocity of 
the post-shock material next to the SNR shell because the 
OH(1720 MHz) masers are only seen in orientation where 
they are transverse to the motion of the expanding shock 
front. Hence they closely match the systemic velocity of the 
SNR. rather than any local intrinsic expansion. We there¬ 
fore expect that molecular gas with similar velocity origi¬ 
nates from the molecular cloud associated with G359.1-0.5. 
In Fig. ^ we show the spectra of 12 CO, HCO + , 13 CO J=l-0 
and CS J=2-l towards maser A. In the velocity range of in¬ 
terest, between —20 and 0 kms -1 , the shape of the 13 CO line 
profile differes from those of the other profiles. Inspection of 
the velocity cubes implies that emission within this range is 
most likely associated with the SNR shell. In contrast, the 
emission between —40 and —20kins -1 and the narrow fea¬ 
tures at +5 and +16 kms -1 appear to be uncorrelated with 
the SNR shell or the H 2 bar. 

Compared with the 13 CO spectra, the 12 CO and HCO + 
spectra, and to a lesser extent the CS spectra, all show an 
absorption feature between —20 and —5kms _1 . In Table [l] 
we list the results of a Gaussian fit to the 13 CO feature, 
but for the other species affected by absorption we give only 
the maximum of the observed line profiles, which thus rep¬ 
resent minimum T m b for these species. The correspondence 
between the peak of the 13 CO emission and the centre of the 
absorption dip in the other three spectra implies that part of 
the molecular cloud is being self absorbed in the molecular 
species with higher optical depth. 

The overall distribution of 12 CO, 13 CO and CS emission 
show differences due to variations in optical depth across 
the region (the spectra of HCO + were not obtained over 
the whole region, as the other molecules, but only towards 
the maser and in four positions around the maser). Maps 
of the emission, integrated between —20 and Okms -1 , are 
shown in Fig. M overlaid with the contours of H 2 emission. 
The true distribution of gas column density is represented 
by the 13 CO because of its low optical depth, in contrast to 
the high optical depths of the other species. We note that 
the 13 CO distribution is dominated by the gas located to the 
west of the H 2 bar. Although the 12 CO and CS emitting 
gas is covering the same region as shown by the 13 CO map, 
these molecules show regions dominated by optically thin 
gas which is not self-absorbed, and is located east from the 
H 2 bar. East of the bar we thus have the gas with lower 
optical depth, while to the west the gas has higher optical 
depth. The 12 CO and CS maps both show an elongated 
feature near and parallel to the H 2 emission and the SNR 
shell. The H 2 bar is located at the edge of this filament of 
molecular gas. 


3.2.2 Physical properties of the molecular cloud 

Self-absorption in molecular clouds can be produced by an 
overlaying colder or sub-thermally excited layer. To esti¬ 
mate the kinetic temperature and density of the gas we have 
modeled the observed values of the 13 CO J=2-l and J=l- 
0 line intensities using a statistical-equilibrium excitation 
code supplied by J. H. Black. The code employs a mean- 
esca pe probability (MEP) approximation for radiative trans¬ 
fer (lansen, van Dishoeck & Black 1994) and calculates line 
intensities given kinetic temperature and density of the gas, 
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Figure 5. 12 CO, HCO+, 13 CO J=l-0 and CS J=2-l spectra 
towards maser A. 
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Figure 6. Maps of 13 CO J=l-0 (upper), 12 CO J=l-0 (middle) 
and CS J=2-l (lower) emission (shown in greyscale with matching 
contours) integrated between —20 and Okms -1 overlaid with the 
contours of H 2 emission. The H 2 contours are: 13, 20, 30, 45, 56 
and 65 xl0 —6 ergs -1 cm -2 sr — T Contour levels (in K km s '' ) 
are: ( 13 CO): 9, 11, 13, 15, 17, 18 and 19; ( 12 CO): 70, 80, 90, 100, 
110, 120 and 130; (CS): 4, 5, 6, 7, 8 and 9. The cross marks the 
location of the OH(1720 MHz) maser A. 


the total column density of the molecule and the line width. 
For the western part of the cloud we find gas temperature of 
8 ± 2 K, gas density of 10 4 cm -3 , 13 CO column density of 
2.0xl0 le cm -2 and 13 CO J=l-0 optical depth of 0.3. This 
low optical depth is consistent with the non-detection of 
C ls O transitions. Also, the density of 10 4 cm -3 is consis¬ 
tent with the molecular distribution being affected by ab¬ 
sorption in this part of the cloud, as the lower-excitation 
temperatures occur at these densities for the HCO + and CS 
lines. Using a 12 CO / 13 CO abundance ratio of 30 ( Gardner 

& Whiteoak 1982) we have the 12 CO column density of 

- =3 - 


O.OxlcF 7 " cm 2 and i2 C0 J=l-0 optical depth of 9. On the 
other hand, the eastern part of the cloud must have density 
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> 10 5 cm -3 because of the CS J=3-2 detection. For such 

not predict OH column density required to produce OH(1720 

densities the corresponding gas temperature derived from 

MHz) masers (Draine, Roberge, & Dalgarno 1983; |Kaufman 

the code is 10 ± 2 K, 13 CO column is 10 lfci cm -2 and 13 CO 

& Neufeld 1996|). Lockett et al. (1999|) and Wardle (1991) 


J=l-0 optical depth is 0.1. This yields a 12 CO column den- 
sity of 3.0xl0 17 cm -2 and optical depth of 3. 

These results show that western region of the cloud may 
be slightly colder than the eastern region. To investigate 
whether this difference in temperature is sufficient to pro¬ 
duce the observed self-absorption, we use the following: 

T abs = f(T bg )e~ Tfa + f(T fg )(l — e~ Tf9 ), (1) 

where 

/(T) = ^[( e w_i)- 1 _( e ^T-l)- 1 ]. (2) 

T abs is the observed temperature at the velocity of maxi¬ 
mum absorption in 12 CO, Tf g is the excitation temperature 
of the foreground part of the cloud, T bg is the excitation 
temperature of the background part of the cloud and r/ 9 
is the optical depth of the foreground cloud. The observed 
T abs ranges from 2 to 3.3 K across the mapped region and 
T crnb = 2.7 K is the temperature of the cosmic microwave 
background. Using T bg = 10 K and Tf g = 3 for the eastern 
part of the cloud we derive the foreground cloud excitation 
temperature of Tf g = 6 ± 1 K. For the western part of the 
cloud T abs depends only on the excitation temperature of 
the foreground part of the cloud because of the high optical 
depth i n that region (rt„ = 91. Fitting the observed TU. we 


derive t ame Tf a of 6 ±1 K. 


The kinetic temperature of the western part of the cloud 
(8 ± 2 K) is in good agreement with the excitation tempera¬ 
ture of 6 ± 1 K derived for the foreground part of the cloud. 
In othe r words, for the western part of the cloud Tkin=. Ze: 


Since t his is the part of the cloud where 12 CO distribution 
is affected by absorption, we conclude that physically cold, 
not sub-thermally excited layer of the cloud is responsible 
for the self-absorption in the observed spectra. 


4 DISCUSSION 

The H 2 1-0/2-1 S(l) line ratio of ~ 20 is consistent with 
shock excitation, as expected from the presence of the 
OH(1720 MHz) maser. Although UV excitation can, in prin¬ 
ciple, produce the 1-0 S(l) line intensity without violating 
the limit on the 2-1 line when molecular gas of densit y 10 5 
- 10 6 cm~ 3 is exposed to a FU V flux of Go > 10 4 ' 5 (Bur¬ 


ton, Hollenbach & Tielens 1990), there is no OB star within 


1 pc, as evident from the 2MASS point source catalogue, 
there are no Hu regions apparent in the continuum map, 
and no significant emission in the Midcourse Space Experi¬ 
ment ( MSX ) bands at 8.3 and 21.3 pm within an aremin 


of 


the 

Thl 


I t 

- 


H» emission must then originate from a shock 


driven by the SNR. The OH(1720 MHz) maser is expected 
to be l ocated within the co oling gas behind a C-type shock 
wave (Lockett et al. 1999), and its location just inside the 
sharp western edge of the H 2 emission supports the idea 
that the expanding shell of G359.5-0.1 is driving a shock 
into a molecular cloud at that location. Shock models do 


http: //www.ipac.caltech.edu/ipac/msx/msx.html 


suggested that the OH column density is enhanced by UV 
photodissocaition of H 2 O. The internal FUV field generated 
by the X-rays from an SNR interior are able to dissociate 
more than 1 per cent of H 2 O, which is e nough to produce the 
required OH abundance (Wardle 1999). How does this FUV 


field affects the H 2 emission? It is about 100 times weaker 
than the standard interstellar field (1.6xl0 -3 ergs -1 cm -2 ), 
which makes it far less than that associated with an H 11 re¬ 
gion or PDR. The corresponding local X-ray energy deposi¬ 
tion per particl e is then H^/n ~ 7 x 10 -29 erg cm 3 s -1 (see 
equation 2 from [ M aloney, Hollenb ach, fe Tielens 1996 ). From 
Fig.6a in Maloney et al. (1990) the intensity in the 2.12 pm 


H 2 1-0 S(l) line of 3.2x 10 -J ergs -1 cm -2 sr -1 which is 0.2 
per cent of observed peak H 2 flux. Thus, the FUV field pro¬ 
duced by X-rays has negligible contribution to the excitation 
of the 2.12 pm H 2 1-0 S(l) line. 

The H 2 bar follows the distribution of extended 


OH(1720 MHz) maser emission (Yusef-Zadeh et al. 1995) 


which is believed to be an evidence of this interaction on the 


global sc ale ( Yuscf-Za dch ct al. 199£ ). Indeed, X-ray obser¬ 
vations ( Bamba et al. 2000| ) show a centre-filled morphology 
of this SNR which, in conjunction with the shell morphol¬ 
ogy in radio band, has been proposed as a characteristic 
of SN Rs interacting with molecular clouds (Rho & Petre 


1995). We found the H 2 emission nea r other masers in thi s 
SNR. aligned well with the SNR shell ( |Lazendic et al. 2001 ). 


The magnetic field along the line of sight towards maser A 
was estimated to be ~ 560pG an d the magnetic field was 
found orient ed along the H 2 bar (Yusef-Zadeh, Roberts & 


Bower 2001). The thickness of the shock is expected to be 


of order 10 17 cm, about an arcsec at a distance of 8.5 kpc. 
We identify the sharp western edge of the bar as being the 
leading, limb brightened edge of a curved shock front. The 
maser velocity is distinct from the mean velocity of about 
—30 km s -1 of the shocked H 2 at this location, but the maser 
is preferentially produced when the shock front is perpen¬ 
dicular to the line of sight. The H 2 velocity gradient could 
be the result of a pre-existing gradient, since a gradient of 
~5kms -1 is found in the 13 CO gas in same direction. 

The results derived from the molecular spectra imply 
the presence of a density gradient across the region, and 
the obvious mechanism to produce such a gradient is the 
passage of an SNR shock wave. The denser gas (~10 5 
10 6 cm -3 ) then corresponds to post-shocked gas and the 
lower density gas (~10 4 cm -3 ) corresponds to ambient pre¬ 
shocked gas, with the H 2 filament being the boundary be¬ 
tween the two regions. A similar case of self absorption 


by pre-shocked gas was also seen in IC 443 (White et al 


1987; van Dishoeck, Jansen, & Phillips 1993). The values ob¬ 


tained for gas densities are consistent with those expected 
in shocks. Because of the presence of the OH(1720 MHz) 
maser, which requires a density ~10 5 cm -3 , the pre-shock 
density must be ~10 4 cm -3 since the typical compression 
in shocks is ~20-30 (for shock speed of 20-30 km s -1 and 
Alfven speed in molecular clouds of 2kms -1 ). In Fig. |we 
present a schematic model of the SNR shock passing through 
the molecular cloud. Dense gas is located mostly on the east¬ 
ern side of the H 2 filament, i.e. behind the shock front, while 
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clouds are found with temperatures low as this (Dickman 


shocked gas 
(molecular hydrogen 
emission) 



1975; Snell 1981). In our case, a likely explanation for this is 
the broad line width of the cloud (15 - 20 kms” 1 ) which en¬ 
ables more efficient cooling. This is supported by the results 
of the MEP cod e. As CO is the dominant coolant at densi¬ 
ties ~10 4 end 3 (Goldsmith & Langer 1978), we can estimate 


the total cooling rate using the total power in the rotational 
transition of 12 CO. For a 12 CO column density of 10 18 cm -2 
and volume density of 10 4 cm -3 , the power emitted for Av = 
20kins -1 and Tki n = 8 K is ~4xl0” 23 ergs” 1 CO” 1 , and for 
Av = 2 kms” 1 and Tki n = 15 K is ~4.7x 10” 23 ergs” 1 CO” 1 . 
Thus the cooling rate from a 8 K cloud with 20kms _1 line 
width is similar to that of a 15 K cloud with a 2 kms” 1 line 
width. 

Despite the positional coincidence and spatial proxim¬ 
ity of the two, we did not find any morphological evidence 
in our data that supports the interaction between the SNR 
and the Snake. Similarly besides the positional coincidence, 
there is no obvious relationship between the H 2 emis sion and 
the Snake. In the view of the recent suggestion by Bickncll 


& Li (2001) that the origin of the Snake is linked to rotat¬ 
ing molecular clouds, we examined the possibility that the 
velocity gradient in H 2 of ±20kms” 1 might originate from 
rotation of a flattened cloud. The 1-0 S(l) emission runs 
along an elliptical depression in the 6-cm emission from the 
SNR where the Snake crosses it (see Fig ^J). This depression 
could be caused by a rotating dense molecular cloud. How¬ 
ever, we derive the mass of the cloud from our 13 CO data 
to be ~ 700 Mq, in which case its escape speed is about 
1.6kms” 1 and its self-gravity could not sustain the rota¬ 
tion. 


5 CONCLUSIONS 


Figure 7. A schematic model of the molecular cloud towards 
the maser A (marked with a cross) is presented to illustrate the 
distribution of different components. The H 2 bar is located just 
behind the thin shock front, which is expanding in a direction 
transverse to the line of sight. 


the lower density gas extends mostly to the west from the 
H 2 emission, i.e. in front of the shock front. 

The gas temperature between the pre-shock and post¬ 
shock gas does not differ much since the gas cools rapidly be¬ 
hind the C-type shock. We did not detect the warm molec¬ 
ular gas, other than in H 2 , found in other SNRs associ¬ 


ated w ith OH( 1 72Q MHz) masers (van Dishoeck, Jansen 
& Phill ps 1993; Frail & Mitchell 1998; Reach fe R.ho 1999| ). 
The reason is probably beam dilution. The warm gas is less 
than 0.6 pc away from the shock front, as inferred from the 
H 2 emission, and our beam of 1-2 pc in size includes the 
volume of gas further away that has been cooled to near¬ 
ambient gas temperature. The post-shock/ambient cloud 
temperature of 10 K is consistent with the excitation tem- 
perature derived for the clouds in the Galactic Centre region 
(Oka et al. 1998). The beam dilution might be also the rea¬ 
son for non-detection of H 2 CO and SiO transitions which 
favour regions with higher temperatures. Although the pre¬ 
shock gas temperature of ~8 K is lo wer than typical molec- 
ular cloud temperatures of > 10 K (Goldsmith 1987), dark 


We carried out near-IR and millimetre wavelength observa¬ 
tions of the region around the OH(1720 MHz) maser located 
at the apparent intersection of the nonthermal filament the 
Snake and the SNR shell G359.1-0.5. A bar of H 2 emis¬ 
sion was found encompassing the OH(1720 MHz) maser and 
aligned with the SNR shell. We suggest that H 2 emission 
originates from the expansion of the SNR blast wave, which 
is evident from the sharp western edge that corresponds to 
the forward shock. This shocked H 2 emission, as inferred 
from the 1-0 and 2-1 S(l) line ratio, supports the notion 
that the OH(1720 MHz) masers association with SNRs are 
produced in molecular shock waves. 

Emission from the molecular species 12 CO, CS, HCO + 
and 13 CO was detected at the maser velocity. The spectra of 
the former three species are affected by broad-line absorp¬ 
tion which originates from a colder layer of molecular gas. 
Optically thin 13 CO was unaffected by the absorption and 
is produced mostly in this colder layer, identified as pre¬ 
shocked gas. The inferred density gradient across the region 
places higher density post-shocked gas west of the H 2 bar. 
The distribution of the pre-shock and post-shock gas is con¬ 
sistent with the passage of the SNR shock and the H 2 bar 
being the boundary between the two regions. 

The warm molecular gas from post-shocked gas was 
not detected in millimetre wavelengths probably because of 
beam dilution. For a better morphological and qualitative 
study of the shocked molecular gas towards maser A we 
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plan to obtain observations at sub-millimetre wavelengths, 
which will have improved spatial resolution. Observations at 
other maser positions in the SNR are in progress to obtain 
a global view of the interaction of the SNR G359.1-0.5 with 
its molecular gas environment. 
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